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The properties of the mass-radius curves of relativistic stellar models constructed from an equation 
of state with a first-order phase transition are examined. It is shown that the slope of the mass- 
radius curve is continuous unless the discontinuity in the density at the phase transition point has 
a certain special value. The curve has a cusp if the discontinuity is larger than this value. The 
curvature of the mass-radius curve becomes singular at the point where the high density phase 
material first appears. This singularity makes the mass-radius curve appear on large scales to have 
a discontinuity in its slope at this point, even though the slope is in fact continuous on microscopic 
scales. Analytical formulae describing the behavior of these curves are found for the simple case of 
models with two-zone uniform-density equations of state. 
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> 

(N 
t> 
O 
(N 
O 
oo 
On 

cr 

u 



X 



I. INTRODUCTION 

The structures of spherically symmetric stellar mod- 
els are usually described in general relativity theory in 
terms of the functions p(r), p(r), and m(r): the total 
energy density, the pressure, and the "mass" contained 
within a sphere of radius r. These functions satisfy Ein- 
stein's equation for a static spherically symmetric space- 
time with fluid source, which may be reduced to the fol- 
lowing pair of ordinary differential equations [jj] : ^] 
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Consider the families of solutions to these equations, each 
of which is determined by a different equation of state 
p = p(p). For each equation of state and for each value 
of the central pressure, there exists a unique solution of 
Eqs. that is non-singular at the center of the star 

Thus for each equation of state there exists a one- 
parameter family of stellar models parameterized by p c 
the central pressure of the star. A large class of equa- 
tions of state have stellar models with finite total radii, 
p(R) = 0, and finite total masses M = m(R). The dis- 
cussion here is limited to these equations of state ||. 
The collection of total masses M(p c ) and radii R(p c ) 
for a given equation of state is called the mass-radius 
curve: [M(p c ), R(p c )]. Each equation of state determines 
a unique mass-radius curve, and conversely it appears 
(although the argument [Q falls short of being a rigorous 
proof) that each mass-radius curve determines a unique 
equation of state. Thus there is hope that the high den- 
sity equation of state of neutron star matter may one 



day be determined by measurements of the macroscopic 
mass-radius curve of these stars. 

This paper is concerned with analyzing the features of 
the mass-radius curve for the case of an equation of state 
with a first-order phase transition. Such a phase transi- 
tion may well be a feature of the equation of state of real 
neutron star matter. Pion condensation Q and/or quark 
deconfinement || might well provide the mechanism that 
drives such a phase transition. This paper does not fo- 
cus on the microphysics of the mechanism that may trig- 
ger such a transition, but rather the consequences that 
such a transition might have on the observable macro- 
scopic equilibrium structures of neutron stars. In partic- 
ular this paper investigates how the properties of such a 
phase transition might be read from the structure of the 
mass-radius curve. 

Consider equations of state that are smooth except at 
one value of the pressure pt where the energy density 
undergoes a simple discontinuity 



lim p{jp) < lim p(jp) = p + 
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FIG. 1. Equation of state with i 
at the pressure p — p t . 
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as illustrated in Fig. |l|. It is convenient to parameter- 
ize the magnitude of the discontinuity in the equation of 



1 



state by the dimensionless quantity A: 

P+ - P- 



A 



P- +Pt 
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FIG. 2. Mass- Radius curves for equations of state with first 
order phase transitions. The different curves correspond to 
different values of the parameter A = (p+ — p_)/(p_ + Pt)- 

Figure ^| illustrates the mass-radius curves for the equa- 
tions of state shown in Fig. ^ (simple polytropes p oc p 2 
with a density discontinuity inserted at the pressure pt). 
The mass scale, M t , and radius scale, Rt, used here are 
the total mass and radius of the stellar model with cen- 
tral pressure pt- The quantity A c that is used to scale 
the discontinuity in the equation of state is defined by 



A, = 



3pt 



2(p-+Pt) 



(5) 



(Note that 1/2 < A c < 3/2.) Figure | illustrates that a 
first-order phase transition makes the mass-radius curve 
bend sharply at the critical point (M t , R t ) where the high 
density phase material first appears in the core of the 
star. 

Figure || makes it appear that mass-radius curves have 
finite discontinuities in their slopes at the point where the 
higher density phase material first enters the star. Fur- 
ther, it appears that the magnitude of the discontinuity 
in the slope is determined by the parameter A that mea- 
sures the magnitude of the phase transition. Thus one 
might hope that an expression can be derived which de- 
termines the properties of the phase transition (e.g. the 
value of A) in terms of some features (e.g. the change 
in slope) of the mass-radius curve in a neighborhood of 
the point (M t ,i? t ). This hope is diminished, however, 
on closer examination of these curves. Figure ^ illus- 
trates the same set of mass-radius curves as shown in 
Fig. ^, however, on a much finer scale. Figure ^ shows 
that the slopes of all of the curves are in fact continuous 
at the point (Mt,Rt), except for the special case with 
A = A c . The mass-radius curves for equations of state 



with A > A c reverse direction at the point (M t ,i? t ), 
however their slopes are continuous there. The curves 
of models with strong first-order phase transitions have 
cusps at the critical point. This microscopic continuity 
of the slope makes it impossible to find a purely local re- 
lationship between the properties of the phase transition 
and the magnitude of the macroscopic bend that occurs 
in the mass-radius curves, as illustrated in Fig. 
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FIG. 3. Mass- Radius curves for equations of state with first 
order phase transitions. This figure represents the same stel- 
lar models depicted in Fig. |^ but on a much finer scale. 

The continuity of the slope of the mass-radius curve 
even in the presence of first-order phase transitions was 
first discovered in the context of Newtonian stellar mod- 
els by Ramsey [[?J and Lighthill ||. The corresponding 
result for relativistic models was demonstrated in anal- 
ogy with the Newtonian analysis by Seidov jjj. A more 
complete and somewhat more rigorous derivation of this 
fact is presented in the Appendix here for arbitrary rela- 
tivistic stellar models. These analyses demonstrate that 
a special value of the magnitude of the phase transition is 
A = A c , with A c defined in Eq. (||). For stronger phase 
transitions, A > A c , the mass-radius curve reverses di- 
rection at the critical point, and probably triggers the on- 
set of instability in the stellar models immediately above 
this point. 

The general analysis of Lighthill and Seidov shows that 
the slope of the mass-radius curve is continuous (in al- 
most all cases) even at the critical stellar model where 
the influence of a phase transition is first felt. This re- 
sult, however, raises more questions than it answers. The 
"typical" mass-radius curves displayed in Figs. ||-^| show 
that the phase transition does have a very profound effect 
on the slope of these curves in a very small neighborhood 
of the critical stellar model. How does the phase transi- 
tion change the curvature of these curves on very small 
scales in a neighborhood of the critical model, while leav- 
ing it relatively unaffected on larger scales? In Section 
U a more detailed analysis of the structure of the mass- 
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radius curve in the neighborhood of a critical point is 
undertaken in an attempt to understand this behavior. 
An analysis is given there of the simple case of equations 
of state having two uniform-density zones: p = p— for 
p < Pt and p = p + for p > p t . Analytical expressions 
are derived for the mass-radius curve for these models 
in a small neighborhood of the critical point (Mt,Rt)- 
These expressions show that the phase transition causes 
the curvature of the mass-radius curve to diverge at this 
point, even though its slope is well defined and contin- 
uous there. This singular part of the curvature causes 
the mass-radius curves in these simple models to bend 
on relatively small scales, much like the more realistic 
ones depicted in Figs. ||-||. Unfortunately, the two-zone 
models are too simple to model accurately the behaviors 
of the mass-radius curves of more realistic equations of 
state. A more complicated analysis is needed, but that 
analysis is deferred to a future investigation. 



II. TWO-ZONE UNIFORM-DENSITY MODELS 

The general solution of Einstein's equation represent- 
ing a static spherical uniform-density star was first found 
by Schwarzschild |10|. Let pi denote the constant density 
of the star. Then the general solution to Eqs. can 
be written: 



m(rj = -^-Pi r + a* j 



p(r) + pt = kfi(r) 



1 + 4tt& 4 / rff{r)dr 



(6) 



(7) 



where at, bi and Ci are arbitrary constants, and fi(r) is 
defined as 



Mr) = 1 



8ir n 2di 
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More complicated stellar models composed of concentric 
uniform-density layers may also be constructed by com- 
bining together the basic solutions given in Eqs. 
These laminated models satisfy Eqs. with pi the 

fluid density within a particular layer. The regularity of 
the global solution is assured by choosing Cj to be the 
inner radius of the i th layer, and the constants a% and bi 
to make p(r) and m(r) continuous at r — C{. 

Now consider the stellar models composed of mate- 
rial having a simple two-zone uniform-density equation 
of state: 



P- 



for p < p t , 
for p > p t . 



(9) 



This is the simplest equation of state having a first order 
phase transition. The family of stellar models associ- 
ated with this equation of state is easily obtained from 
Eqs. (§)-©. 



For stars with small central pressures p c < pt, the 
solutions are the standard interior Schwarzschild mod- 
els jioj . These may be obtained from the general expres- 
sions above by setting ejj = c, = and bi = p c + p_ . The 
total masses and radii of these models are determined by 
solving Eqs. (||)-(@) for the points where p(R) = and 
M = m(R). As functions of the central pressure these 
solutions are: 



R 2 (Pc) 



- Pc 



3p c 



M{p c 



-in 



P-R\p c ). 



(10) 



(11) 



Thus the models with small central pressures have the 
well-known cubic mass-radius curve. The critical model 
in this family is the one having p c — p t with total mass 
M t — M(p t ) and total radius R t — R(pt)- Expressions 
for these critical values are given by Eqs. (|lTj|)-(pj]) with 
Pc =Pt- 

The stars with large central pressures, p c > pt, have 
two concentric layers. The inner layer is composed of 
high density material with p = p + , and the outer layer of 
lower density material with p — p_ . The structure of the 
inner core is determined from Eqs. (||)-(@) with pi = p+, 
a i = c i — 0, and bi = p c + p + . The radius, r t , of this 
inner core is determined by solving Eq. ([?]) for the point 
where pt = p(ft)- As a function of the central pressure, 
p ci this core radius is 



(P+ +Pc){p++3pt) 
(P+ +Pt)(p+ + 3p c ) 



(12) 



The outer envelopes of these models are determined again 
by Eqs. (|)-(B) with Pi = p_ , Oi = 4ir(p+ - P-)rf/3, 
bi = (p- +pt)y/l — 87T/9+rf/3, and Ci — r t . The quadra- 
ture indicated in Eq. (|7) can be expressed in terms of 
standard elliptic integral functions for this case, but that 
representation does not offer any particular insight for 
our purposes here. 

The total masses and total radii of the two-zone 
uniform-density models are found by solving Eqs. 
(H) for the points where p(R) = and M = m(R). 
These equations can not be solved analytically even for 
these simple two-zone models. However, the nature of 
the solutions near the critical model can be studied by 
means of power series expansions. The small parameter 
s = (r t /Rt) 2 (which vanishes as p c J, p t ) can be used to 
expand the various quantities (i.e., bi and fi) that ap- 
pear in Eq. ([?]). The integration that appears in Eq. (Q) 
is performed term by term and the resulting equation is 
solved for the total radius of the star p(R) = 0. The 
resulting series expansion for R is used to evaluate the 
total mass M = m(R) using Eq. (|^) to the same order of 
approximation: 
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Rt 



= 1 + ns + r 3/2 s 3 / 2 + r 2 s 2 + 0(s 5 / 2 ), (13) 



5/2N 



M(s) 
M, 



= 1 + mis + m 3/2 s 3/2 + m 2 s z + 0(s o/z ), (14) 



5/2n 



where the expansion coefficients r*i, mi etc. are given by 
A c - A 

n = ! (15) 

(16) 



^3/2 



8A3 ' 
A(3 - 6Ag - 8A^) 
8A4(3-2A C ) : 



^2 = 



9A(4A C -A)(4A 2 -1) r 2 
128A6(3 - 2A C ) ~2 
mi = 3r 1: 

_ A(9 - 18Ag - 8A 4 C ) 
1713/2 ~ 8A4(3-2A C ) ' 
m 2 = 3(r 2 + r 2 ), 



(17) 
(18) 
(19) 
(20) 



and where A and A c are defined in Eqs. (||) and (||) 
respectively. These series expansions give a reasonably 
good approximation of the mass-radius curve near the 
critical model as illustrate in Fig. |[ The series agree with 
the exact (numerically determined) mass-radius curves to 
within 1% for models whose masses and radii differ from 
the critical values by up to about 10% in the worst case 
examined. (The series converge most poorly for A = A c 
among those cases examined.) This level of accuracy is 
good enough to account for the interesting features of 
the mass-radius curve near the critical point. Figure |] 
illustrates that the series exhibit the same apparent dis- 
continuity in slope on large scales as the exact curves. 




R/R t 

FIG. 4. Mass-Radius curve for the two-zone uni- 
form-density equation of state. The series expansion for the 
curve (dashed line) is compared to the exact (solid line) for a 
model with A = 2A c /3 and Pt/p- = 0.5. 

The tangent vector to the mass-radius curves is de- 
termined by differentiating the expressions for x(s) = 
R{s)/R t and y(s) = M(s)/M t in Eqs. P)-©: 



dx/ds 
dy/ds 



ri + §r 3 / 2 s 1/2 



mi 



2r 2 s 



§m 3 / 2 s 1/2 + 2m 2 s 



(21) 



This expression illustrates that the slope of the mass 
radius curve is continuous even at the critical model. 
This follows from the fact that the tangent vector com- 
puted from just above the critical model using Eq. ( ^j~Lj ) is 
ri(dx/dp c )^ 1 times that computed from just below using 
Eq. ([ll]) . Just above the critical point the tangent vector 
is proportional to the quantity r\ defined in Eq. ( |l5| ) . r% 
is positive for weak phase transitions A < A c but neg- 
ative for strong transitions A > A c . Thus the tangent 
vector and hence the mass-radius curve itself reverses di- 
rection when a strong phase transition occurs. This con- 
firms the general continuity analysis of Ramsey, Lighthill, 
and Seidov for the simple case of the two-zone uniform- 
density models. 

While the general analysis of the structure of the mass- 
radius curve fails for the case of a phase transition with 
A = A c , the analysis here for the simple two-zone 
uniform-density models succeeds even in this case. The 
terms proportional to s in Eqs. (|D|)-(^|) vanish when 
A = A c , and therefore s is not a good affine parame- 
ter for the mass-radius curve at s = in this case. In- 
stead, the appropriate parameter is A = s 3 / 2 . In this 
case the tangent vector evaluated at the critical model 
is (dx I dX, dy / d\) — (?"3/ 2 , m 3 / 2 ). This vector is not pro- 
portional to the tangent vector just below the critical 
model, and so the slope of the mass-radius curve is not 
continuous in this case. 

In order to discuss the magnitude of the change in slope 
that occurs at the critical model it is necessary to adopt 
a metric structure for the space of masses and radii. This 
makes it possible to define the inner products of tangent 
vectors (and so define angles) and also more generally to 
discuss the curvatures of these curves. There is no canon- 
ical choice for this metric, and therefore no absolute in- 
trinsic meaning can be given to angles or curvatures that 
are computed. Nevertheless these quantities are useful 
tools for understanding the features of the mass-radius 
curves seen in Figs. g-Q. Thus the metric used to display 
those figures is adopted: the flat metric with Cartesian 
coordinates x = R/Rt and y = M/M t . 

Return now to the kink that occurs in the mass-radius 
curve for the case of a phase transition with A = A c . The 
angle between the slopes above and below the transition 
point can be determined (using the metric defined above) 
from the inner product between the tangent vectors. The 
resulting angle 9 depends solely and monotonically on the 
ratio pt/p~: 



cos 6 = 



r 3 /2 + 3m 3/2 
10(i 



3/2 



'3/2 



(22) 



This 8 varies from about 4.4° for pt/p- — to about 
170.4° for Ptj P- — 00. This formula is a simple example 
of the kind of relationship that one had hoped to find 
relating the parameters of the phase transition and the 
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macroscopic structure of the mass-radius curve. In this 
special case (phase transitions with A = A c ) the magni- 
tude of the kink in the mass-radius curve determines the 
ratio pt I p- ■ Unfortunately, this formula is not universal 
even for phase transitions with A = A c . The magnitude 
of the kink displayed in Fig. || does not satisfy this equa- 
tion for example. The general form of this relationship 
must depend on other features of the equation of state 
(e.g. dp/ dp at the transition point) that are not present 
in the simple two-zone uniform-density models. 

The curvature of any of the mass-radius curves can 
be evaluated by differentiating the unit tangent vector 
along the trajectory of the curve. The resulting acceler- 
ation is equal to the inverse of the radius of curvature of 
the curve. For a general curve in a flat two-dimensional 
space, this acceleration is given by 



d 2 x dy 
ds 2 ds 



d 2 y dx\ 
ds 2 ds J 



dx\' 
~ds) 



-3/2 



(23) 



This expression is invariant under changes in the param- 
eterization of the curve, but not on the assumed metric of 
the mass-radius space. It is straightforward to evaluate 
this expression using the series expansions, Eqs. (pX 
(fj), for the curve: 



3(3r 3/ 2 - m 3/2 ) 
40 
3 



5V10 



lOr? 
0(s^ 2 ). 
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1TJ2 



9(^3/2 + 3to 3/2 ) 
20ri 



(24) 



The first term in Eq. ([24]) is proportional to A and there- 
fore vanishes when there is no phase transition. The sec- 
ond term is a pure number that is independent of the 
parameters of the phase transition. The second term is 
therefore the curvature of the mass-radius curve just be- 
low the phase transition point. The first term includes 
a part proportional to 1/s 1 / 2 which diverges at the criti- 
cal model. This infinite curvature causes the mass-radius 
curve to bend sharply in the neighborhood of the critical 
point. 

The analysis here shows that the mass-radius curves 
of stellar models with first-order phase transitions have 
infinite accelerations at the critical model where the high 
density phase first appears. This acceleration causes 
these curves to bend sharply, appearing on large scales 
to have a discontinuous slope at this point. Analytical 
formulae, Eqs. (|l^)-(|l4|), have been derived that describe 
quantitatively the structures of these curves for models 
with simple two-zone uniform-density equations of state 
as illustrated in Fig. ||. These formulae also account 
in a qualitative way for the behavior of the mass-radius 
curves of more realistic equations of state, as illustrated 
in Figs. ||-^. The quantitative description of the more 
realistic mass-radius curves (as would be needed to ana- 
lyze the measured masses and radii of real neutron stars) 
requires the derivation of the analogs of Eqs. (ft3l)-dlT 



for a general equation of state. This generalization is not 
a simple extension of the analysis presented here, and is 
deferred to a future investigation. 
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APPENDIX: CONTINUITY OF THE SLOPE OF 
THE MASS-RADIUS CURVE 

This appendix presents the argument that the slope 
of the mass-radius curve is continuous even at the point 
where the phase transition first sets in (except for the 
case A = A c ). The discussion here is more complete and 
somewhat more rigorous than Seidov's original ||. 

This argument can be made a little simpler by intro- 
ducing a somewhat unusual representation of the equa- 
tions of stellar structure. The structure of a spherical 
star in general relativity is usually expressed in terms 
of the functions m(r, p c ), p{r,p c ) and p(r,p c ) satisfy- 
ing Eqs. (Q)-(^). The analysis here is concerned with 
how these functions behave as p c varies near the point 
Pt where a phase transition occurs. Unfortunately the 
functions m(r,p c ), p(r,p c ), and especially p(r,p c ) are not 
smooth when a phase transition is present. The den- 
sity function p(r,p c ) in particular is discontinuous at the 
boundary of the inner core of high density phase material, 
and the position of this discontinuity changes as p c is var- 
ied. Thus, approximate expressions for these functions in 
terms of power series expansions near the critical point 
(as needed in Seidov's analysis) are somewhat awkward. 
This difficulty is simplified by considering the structure 
of the star in terms of the "inverses" of these functions 
Q: e.g. m(p,p c ) and r(p,p c ). Since the pressure is a 
monotonic function of the radius in these models, this 
inversion is always possible. These functions are smooth 
in their dependence on p for fixed p c , and so it is more 
straightforward to approximate them with power series 
expansions. 

It is also useful to introduce a slightly different set of 
basic variables to describe the structures of stars instead 
of the usual p, m, and r. It is preferable to use the 
thermodynamic enthalpy function 



Hp) 



dp 
o P + P 1 



(Al) 
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in place of the pressure as the independent variable in 
this representation of the problem, because it makes the 
differential equations non-singular at the surface of the 
star. Similarly, it is somewhat preferable to use the func- 
tions u = r 2 and v = m/r as dependent variables because 
they are smoother functions of h near the centers of the 
stars. The straightforward translations of the standard 
structure equations (|l])-(||) into this new set of variables 
gives 



du 
dh 



2u(l-2«) 

= U {u, v, h), 



4:Ttup(h) + v 



(A2) 



— = -(l-2v) \\,, =V{u,v,h). (A3) 
an VKUpyti) + v 

In these Eqs. (|A2|)-([A^) the functions p(h) and p(h) 
are determine d fr om the chosen equation of state p = 
p{p) and Eq. (Al). They are therefore explicitly known 
functions once a particular equation of state has been 
selected. This version of the equations has several nice 
features. First, the use of h as the independent variable 
makes the domain where the solution is defined, [0, h c ] 
where h c is the value of h at the center of the star, known 
before the solution is found rather than after. Second, the 
total radius of the star is determined simply by evaluating 
the function u at the surface of the star h = 0, instead 
of solving the usual surface equation p(R) = 0. Third, 
the use of h as independent variable makes the equations 
non-singular at the surface of the star. And fourth, the 
use of u and v as dependent variables make the solutions 
near h — h c smoother than the usual functions m and r. 

Consider the one-parameter family of solutions to these 
equations constructed from a single equation of state: 
u(h,X) and v(h, A), where A is the parameter that dis- 
tinguishes the individual members of the family. Each 
member of this family satisfies the usual boundary con- 
ditions, both at the center of the star h = h c , 



u[h c (X),X]=v[h c (\),X]=0, 
and at the surface of the star h = 0, 
u(0,A) = R 2 (X), 



v(0,X) 



M(X) 
R(X)> 



(A4) 



(A5) 



(A6) 



where M(A) is the total mass and R(X) is the total ra- 
dius of the model with parameter A. The choice of the 
particular parameterization is arbitrary; however, it is 
convenient to insist that each member of the family have 
a unique central pressure p c and a unique central value 
of the enthalpy h c . Thus, either of these quantities could 
be used as the parameter A. 

Near the centers of these stars, the solutions to the 
structure equations can be given analytically as power 



series expansions. When the equation of state is smooth 
[i.e., when p(h) and p(h) are smooth functions] then 
u(h, A) and v(h, A) have the expansions: 



u(h, A) 



v(h, A) 



3(/t e - h) 
2n(p c + 3p c ) 



Zpcjhc - h) 
Pc + 3p c 



+ 0{h c -hf 



+ 0(h c -h) 2 



(A7) 



(A8) 



The right sides of Eqs. (|A7|)-(|Aq) depend on A implic- 
itly. The central value of h depends on which member 
of the one-parameter family is being considered, thus 
h c = h c (X). The choice of parameterization is arbitrary 
however. Thus h c (X) is an arbitrarily monotonic func- 
tion. The quantities p c and p c also depend on A in the 
obvious ways: p c = p[h c (X)], etc. 

Next consider the situation where the equation of state 
is smooth, except at a certain phase transition point h = 
h t . Assume that the density has a finite discontinuity at 
this point: 



9- 



lim p(h) < lim p(h) 

h1h t h[h t 



P+- 



(A9) 



This is simply the restatement of Eq. (|3|) in terms of h in- 
stead of p. The pressure f unc tion p(h) is C° at this point 
as a consequence of Eq. (|Al| ) , but it has a finite discon- 
tinuity in its first derivative there. One of the important 
facts that is needed in this analysis is the continuity of 
the functions u and v. The functions u(h, A) and v(h, A) 
are C° functions of h for fixed A and C° functions of A 
for fixed values of h. These continuity conditions are rea- 
sonably easy to establish. First consider the continuity 
of u(h, A) and v(h, A) as functions of h for fixed A. The 
right sides of Eqs. (A2 )— (|A3| ) are smooth functions of u 
and v (for u > and v > 0) and are continuous func- 
tions of h e xcep t when h = h t . When h = ht the right 
side of E q. (A3) has a finite discontinuity as described in 
Eq. ( |A9| ) . If u and v were disc onti n uous for some value of 
h, then the left sides of Eqs. (A2)-(A3) would be singu- 



lar there. But the right sides are finite for h < h c (X), so 
u(h, A) and v(h, A) must be continuous for all h < h c (X) 
for fixed A. 

Next consider the continuity of u[h, A) and v(h, A) as 
functions of A for fixed h. Assume that the parameter 
A is chosen so that h c (X) is smooth and monotonically 
increasing. Let At denote the critical value of the pa- 
rameter for which h c (Xt) = h t . For A < At the expan- 



sions in Eqs. (A7)-(A8) show that u(h,X) and v(h, A) 
are continuous in A at least in a small neighborhood of 
the center of the star where h = h c . The differential 
Eqs. ( A2)-([A3|) are non-singular outside of this neigh- 
borhood. The standard theorems insure that the 
solutions to such non-singular equations depend contin- 
uously on their boundary values. These boundary values 
as determined by Eqs. ( A7)-(|Ag) can be applied a small 
distance away from the singular point h — h c . Thus, 
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u(h,X) and v(h, A) are continuous functions of A for fixed 
h, at least for A < A f . When A > A t a similar argument 
insures the continuity of u(h, A) and v(h, A) in the cores 
of the stars where h > h t . Thus u(h t ,X) and v(h tl X) 
are continuous functions of A for A > A*. These func- 
tions can now be considered as the boundary values for 
u(h, A) and v(h, A) in the domain h < ht- In this domain 
the standard theorems again apply, so the continuity of 
the boundary values [i.e., u(h t , A) and v(h t ,X)} guarantee 
the continuity of u(h, A) and v(h, A) as functions of A for 
fixed h. The only troublesome point is at A = A t . 

Consider the stellar models with A just above the crit- 
ical point. These models consist of a very small central 
core of material of the higher-density phase, p> p+, and 
the vast majority of the material in the lower-density 
phase. In the limit A J. At the size and mass of this cen- 
tral core of material goes to zero. This limit c an be see n 
analytically in the expansions given in Eqs. (|A7|)-([A8|). 
In this limit what remains is a star composed entirely 
of matter in the lower-density phase, except for the sin- 
gle point at the center of the star. At this single central 
point the material remains in the higher-density phase. 
But, the matter at this single point does not effect the 
struct ure o f the star at all. The solutions to the structure 
Eqs. (|A2[)-(|A3|) are not changed if the equation of state 
is changed only at a single value of h. Thus the func- 
tion ui(h) = liniAXA t u(h,X) is identical to the function 
that describes a stellar model consisting entirely of lower 
density material with h c = h t : u-\ = \mi\^\ t u(h, A). A 
similar argument applies to v(h, A). Thus the functions 
u(h, A) and v(h, A) are continuous functions of A even for 
A = A t . 

In order to understand the structure of the mass-radius 
curve in stars having a first-order phase transition, the 
structure of stars having a very small central core of the 
high density phase material must be analyzed in some 
detail. The central cores of such mode ls are described by 
the series solutions given in Eqs. (A7)-(A8): 



u+(h,\) = 



3{h c - h) 
2n(p c + 3p c ) 



+ 0{h c -hf 



v + (h,X)= 2pc{k c k) +0(h c -hf, 
Pc + 3p c 



(A10) 



(All) 



for ht < h < h c (X). The outer envelopes of these stars 
are composed of material from the lower density phase. 
In the stars of interest here — those with only very small 
cores of high density material — the structure of the outer 
envelope is nearly identical to the structure of a star com- 
posed entirely of low density phase material. Thus the 
inner region of this outer envelope may be approximated 
as 



U-(h, A) = 



3(/i c - h) 



Mp- 

+0{h c 



8u{h) (h c — h t ) 



(A12) 



v. (h, A) = 2p - {h \ k) + Sv(h) (h c - h t ) 
p- + 3p c 

+0(h c -h) 2 , 



(A13) 



where 5u and 5v are solutions to the linearized struc- 
ture equations. Quite generally, these linearized struc- 
ture equations have the form 



dSu dU r dU r 

IT = ~x~ du + ~7T dv > 
ah ou av 



ddv 
~dh 



dV 
du 



5u 



dV 
dv 



5v, 



(A14) 



(A15) 



where U(u.v, h) a nd V(u, v, h) are the functions defined 
in Eqs. (|A2|)-(A3). For our purposes here it is sufficient 
to evaluate these functions using the first orde r term s in 
the expansions for u and v given in Eqs. (A7)-(A8). In 
this case the functions of interest to us have the forms: 



— = 2n(p_+ Pc ) — 
_ 27r(/9_ +p c )p- 



p- + 3p c 



8V „ . ,dU 
— = 2tt(p_ +p c )^r- 
ov ov 

S(p-+ Pc ) 1 



+ 0{h c ~ h)°, (A16) 



2(p_ + 3p c ) h c — h 



+ 0(h c -h) . 



(A17) 



The resulting form of Eqs. (A14)-( A.15) can be integrated 
analytically. The general solution is 



Su(h) = A + 



D 



6v{h) 



Airp- 



A ■ 



{h c -h)V* [ 



2?r(p_ +p c ) 

{h c - hy/ 2 



B. 



(A18) 



(A19) 



where A and B are arbitrary constants. The values of 
these constants are determined by demanding continuity 
at h = h t of the functions describing the inner core, u+ 
and v + , with the functions describing the outer envelope, 
it_ and These continuity conditions are satisfied for 
the following values of A and B: 



A = 



B = 



Hpc-p-) 



2i:(p c + 3p c )(p- + 3p c ) ' 



3{ Pc ~ p^){h c -h t ) 1 / 2 
n(p c + 3p c ){p- + 3p c ) ' 



(A20) 



(A21) 



The complete expressions then for the inner portions of 
the structure functions in the low-density envelope of the 



star are 
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U-(h, A) 



3{h c -h) 



9(p c - P-)(h c - h t ) 



2-rt(p- + 3p c ) 2n(p c + 3p c )(p- + 3p c ) 

0{h c -hf, (A22) 



2 

l -3 



h r — h t 



1/2 



v-(h,X) 



2p_ (h c - h) _ 6(p c - p-){h c - h t ) 
p-+3p c {p c + 3p c )(p- + 3p c ) 



P- - (P- +Pc) 



+0(h c - hf 



h c - h t 



1/2 



(A23) 



The match of u_ to the inner-core function u + is C 1 at 
h = ht, as required by Eq. (A2). The match of u_ to v + 
is C° at h = ht- The slope of V- differs from that of v+ 
at h = h t by the amount required by Eq. (A3). 

Now consider the region in these models where \h c — 
h t \ <C \h c — h \ <CL The approximate expressions given 



m Eqs. (|A2g)-(|A2f) are valid for these models in this 
region. In addition the terms in these expressions pro- 
portional to (h c — ht) 1 ^ 2 can be neglected in this region 
as well. Next, evaluate the derivatives Su^ — du/dX, etc. 
in this region for these models having a very small core 
of high-density phase material: 



8u±(h, A) 



d u_ 
~~dX~ 



3[-2p c + 3(p_+p c )] dh c 



2tt( P - + 
+0(h c - h) + O 



3pc)(p c 
h c - h t 



3p c ) dX 

1/2 



(A24) 



5v±(h, A) 



dv- _ 2p_[-2p c + 3(p_ + Pc )] dh c 



OX 



-0{h c 



(P- 
h)- 



3pc)(p c + 3p c ) dX 



O 



h r — h 



1/2 



(A25) 



These expressions can now be compared with those that 
pertain to stars having no material at all in the high- 
density phase. Thus, define <5wj = du/dX, etc. for the 
models with no high d ensi ty pha se material at all using 
the expansions in Eqs. ( A7)-(|Aq) that are valid through- 
out the inner regions of these models: 



Sui(h, A) = ^ 



dh r 



2tt(p- + ip c ) dX 



+ 0(h c 



h), 

(A26) 



5vf(h, A) = 



2p_ 



dh r 



dv 



dX p- + 3p c dX 



OQi c -h). (A27) 



These expressions illustrate that the derivatives 5u and 
6v are not continuous functions of A near the critical 
model with A = A*. However, the discontinuity is of a 
very special type. These expressions for (Sui,Svi), are 
related to those for (Su^, Sv^) in the following simple way, 



5u l 



= 2 



A c - A 
2A r + A 



5uf 



(A28) 



where A and A c are d efined in Eqs. (^) and (||) respec- 
tively. Equation ( |A2§| ) is exact in the limit that A — > At 
from above and below respectively, and when the func- 
tions are evaluated at the center of the star h = h t . 
The derivatives Su and Sv satisfy the linear differen- 



tial Eqs. (A14)-(A15). Further, the continuity of u and 
v as functions of A at the point A = A t , implies that 
both (<5itj, (5vf ) and (6u±,6vi) satisfy the same differen- 
tial equation at A = A t . (This fact is the reason that it 
was necessary to establish the continuity of u and v in 
some detail above.) Thus, it follows that the functions 
(<5iij, <5i>l) are proportional to (du±,6vi) throughout the 
critical model with A = At since they are proportional to 
one another in a neighborhood of h = h t . 

At the surface of the star, h = 0, the derivatives 5u 
and Sv are related to the tota l ma ss and radius of the 
star as a consequence of Eqs. ( A5)-(|A6|). In particular 
these functions must satisfy: 



Su(0,X) = 2R{X) 



dR(X) 
dX ' 



6v{0, A) 



1 dM(X) M(X) dR(X) 



R(X) dX R 2 (X) dX 



(A29) 



(A30) 



The functions (Su^, Sv^) (evaluated for models just above 
the critical one) are proportional to (5u^, 5v-\) (evaluated 
for models just below the critical one) throughout the 
star. Thus, the surface values of these functions are pro- 
portional as well. This implies in particular that the 
tangent vectors to the mass radius curve evaluated above 
and below the critical model are related by: 



(A31) 



This expression has several interesting consequences. 
First, it shows that the tangent vector to the mass- 
radius curve is discontinuous at the critical model when- 
ever one parameterizes the curve in a way that makes 
h c (X) smooth. Second, this expression shows that the 
mass-radius curve in fact reverses direction at the criti- 
cal model if the phase transition is su fficie ntly severe so 
that A > A c . Third and finally, Eq. ( |A3l| ) implies that 
the slope of the mass-radius curve, dM/dR, is continuous 
even at the critical model: 




dM _ dM l (dR^ 1 _ dM 1 (dR^ ' 
~d~R ~ ~dX~ \~dX 



dX V dX 



(A32) 



Continuity of the slope pertains even if the curve has a 
cusp and reverses direction at the critic al model, unless 
A = A c . In this special case Eq. (A31) merely implies 
that A [chosen so that h c (X) is smooth] is not a good 
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affine parameter for the mass-radius curve. A higher or- 
der analysis is needed to understand the differentiability 
of the curve in this special case. 
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